# -*- coding: mbcs -*-
# Do not delete the following import lines
from abaqus import *
from abaqusConstants import *
import __main__

def deviation(partname, k1, k2, k3, k4 ,k5, k6, k7, k8, k9, k10, k11, k12, k13, k14, k15, k16, k17, k18, k19, k20, k21, k22, k23, k24, k25, k26, k27, k28, k29, k30, k31, k32, k33, k34, k35, k36, k37, k38, k39, k40, k41, k42):
    import section
    import regionToolset
    import displayGroupMdbToolset as dgm
    import part
    import material
    import assembly
    import step
    import interaction
    import load
    import mesh
    import optimization
    import job
    import sketch
    import visualization
    import xyPlot
    import displayGroupOdbToolset as dgo
    import connectorBehavior
    s = mdb.models['Model-1'].ConstrainedSketch(name='__profile__', 
        sheetSize=500.0)
    g, v, d, c = s.geometry, s.vertices, s.dimensions, s.constraints
    s.setPrimaryObject(option=STANDALONE)
    s.ConstructionLine(point1=(0.0, -250.0), point2=(0.0, 250.0))
    s.FixedConstraint(entity=g[2])
    s.ArcByCenterEnds(center=(0.0, 0.0), point1=(0.0, 120.0), point2=(107.5, 45.0), 
        direction=CLOCKWISE)
    s.CoincidentConstraint(entity1=v[2], entity2=g[2], addUndoState=False)
    s.CoincidentConstraint(entity1=v[0], entity2=g[2], addUndoState=False)
    s.Line(point1=(110.692918596157, 46.3365705751356), point2=(130.089622557846, 
        0.0))
    s.TangentConstraint(entity1=g[3], entity2=g[4], addUndoState=False)
    s.ArcByCenterEnds(center=(140.0, 0.0), point1=(130.089622557846, 0.0), point2=(
        155.0, 0.0), direction=COUNTERCLOCKWISE)
    s.ConstructionLine(point1=(149.910377442154, 0.0), point2=(130.089622557846, 
        0.0))
    s.HorizontalConstraint(entity=g[6], addUndoState=False)
    s.CoincidentConstraint(entity1=v[4], entity2=g[6], addUndoState=False)
    s.CoincidentConstraint(entity1=v[3], entity2=g[6], addUndoState=False)
    s.HorizontalConstraint(entity=g[6])
    s.FixedConstraint(entity=g[6])
    s.TangentConstraint(entity1=g[4], entity2=g[5])
    s.RadialDimension(curve=g[3], textPoint=(47.87109375, 53.1080627441406), 
        radius=120.0)
    s.RadialDimension(curve=g[5], textPoint=(114.669555664063, -30.4274921417236), 
        radius=10.7436438404634)
    s.AngularDimension(line1=g[4], line2=g[6], textPoint=(99.1864624023438, 
        2.7215404510498), value=67.2855876468284)
    s.DistanceDimension(entity1=v[5], entity2=g[2], textPoint=(27.5218811035156, 
        -59.5986328125), value=140.0)
    s=mdb.models['Model-1'].sketches['__profile__']
    s.Parameter(name='R', path='dimensions[0]', expression=k1)
    s.Parameter(name='r', path='dimensions[1]', expression=k2, 
        previousParameter='R')
    s.Parameter(name='D', path='dimensions[3]', expression=k3, 
        previousParameter='r')
    s.Parameter(name='a', path='dimensions[2]', expression=k4, 
        previousParameter='D')
    s.delete(objectList=(g[6], ))
    s.undo()
    s.sketchOptions.setValues(constructionGeometry=ON)
    s.assignCenterline(line=g[2])
    p = mdb.models['Model-1'].Part(name=partname, dimensionality=THREE_D, 
        type=DEFORMABLE_BODY)
    p = mdb.models['Model-1'].parts[partname]
    p.BaseShellRevolve(sketch=s, angle=360.0, flipRevolveDirection=OFF)
    s.unsetPrimaryObject()
    p = mdb.models['Model-1'].parts[partname]
    session.viewports['Viewport: 1'].setValues(displayedObject=p)
    del mdb.models['Model-1'].sketches['__profile__']
    p = mdb.models['Model-1'].parts[partname]
    p.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=42.0)
    p = mdb.models['Model-1'].parts[partname]
    p.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=84.0)
    p = mdb.models['Model-1'].parts[partname]
    p.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=126.0)
    session.viewports['Viewport: 1'].view.setValues(nearPlane=766.851, 
        farPlane=1398.39, width=472.177, height=216.022, viewOffsetX=-90.9582, 
        viewOffsetY=-19.0706)
    p = mdb.models['Model-1'].parts[partname]
    p.DatumPlaneByPrincipalPlane(principalPlane=XYPLANE, offset=0.0)
    p = mdb.models['Model-1'].parts[partname]
    p.DatumPlaneByPrincipalPlane(principalPlane=YZPLANE, offset=0.0)
    p = mdb.models['Model-1'].parts[partname]
    d1 = p.datums
    p.DatumPlaneByRotation(plane=d1[5], axis=d1[1], angle=45.0)
    p = mdb.models['Model-1'].parts[partname]
    d2 = p.datums
    p.DatumPlaneByRotation(plane=d2[5], axis=d2[1], angle=135.0)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    pickedFaces = f.getSequenceFromMask(mask=('[#7 ]', ), )
    d1 = p.datums
    p.PartitionFaceByDatumPlane(datumPlane=d1[2], faces=pickedFaces)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    pickedFaces = f.getSequenceFromMask(mask=('[#f ]', ), )
    d2 = p.datums
    p.PartitionFaceByDatumPlane(datumPlane=d2[3], faces=pickedFaces)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    pickedFaces = f.getSequenceFromMask(mask=('[#1f ]', ), )
    d1 = p.datums
    p.PartitionFaceByDatumPlane(datumPlane=d1[4], faces=pickedFaces)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    pickedFaces = f.getSequenceFromMask(mask=('[#3f ]', ), )
    d2 = p.datums
    p.PartitionFaceByDatumPlane(datumPlane=d2[5], faces=pickedFaces)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    pickedFaces = f.getSequenceFromMask(mask=('[#fff ]', ), )
    d1 = p.datums
    p.PartitionFaceByDatumPlane(datumPlane=d1[8], faces=pickedFaces)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    pickedFaces = f.getSequenceFromMask(mask=('[#ffffff ]', ), )
    d2 = p.datums
    p.PartitionFaceByDatumPlane(datumPlane=d2[6], faces=pickedFaces)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    pickedFaces = f.getSequenceFromMask(mask=('[#ffffffff #f ]', ), )
    d1 = p.datums
    p.PartitionFaceByDatumPlane(datumPlane=d1[7], faces=pickedFaces)
    session.viewports['Viewport: 1'].partDisplay.setValues(sectionAssignments=ON, 
        engineeringFeatures=ON)
    session.viewports['Viewport: 1'].partDisplay.geometryOptions.setValues(
        referenceRepresentation=OFF)
    mdb.models['Model-1'].Material(name='Material-1')
    mdb.models['Model-1'].materials['Material-1'].Density(table=((k5, ), ))
    mdb.models['Model-1'].materials['Material-1'].Elastic(table=((k6, k7), 
        ))
    mdb.models['Model-1'].materials['Material-1'].Plastic(table=((k8, 0.0), 
        ))
    mdb.models['Model-1'].ExpressionField(name='AnalyticalField-1', localCsys=None, 
        description='', 
        expression=k11 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-1', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-1', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-2', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-1'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-3', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-2'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-4', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-3'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-5', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-4'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-6', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-5'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-7', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-6'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-8', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-7'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-9', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-8'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-10', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-9'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-11', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-10'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-12', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-11'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-13', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-12'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-14', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-13'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-15', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-14'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-16', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-15'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-17', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-16'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-18', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-17'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-19', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-18'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-20', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-19'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-21', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-20'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-22', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-21'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-23', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-22'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-24', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-23'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-25', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-24'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-26', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-25'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-27', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-26'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-28', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-27'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-29', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-28'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-30', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-29'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-31', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-30'])
    mdb.models['Model-1'].AnalyticalField(name='AnalyticalField-32', 
        objectToCopy=mdb.models['Model-1'].analyticalFields['AnalyticalField-31'])
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-2', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-2', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-3', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-3', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-4', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-4', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-5', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-5', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-6', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-6', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-7', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-7', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-8', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-8', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-9', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-9', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-10', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-10', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-11', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-11', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-12', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-12', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-13', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-13', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-14', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-14', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-15', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-15', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-16', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-16', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-17', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-17', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-18', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-18', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-19', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-19', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-20', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-20', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-21', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-21', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-22', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-22', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-23', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-23', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-24', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-24', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-25', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-25', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-26', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-26', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-27', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-27', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-28', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-28', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-29', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-29', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-30', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-30', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-31', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-31', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].HomogeneousShellSection(name='Section-32', 
        preIntegrate=OFF, material='Material-1', thickness=0.0, 
        thicknessType=ANALYTICAL_FIELD, thicknessField='AnalyticalField-32', 
        idealization=NO_IDEALIZATION, poissonDefinition=DEFAULT, 
        thicknessModulus=None, temperature=GRADIENT, useDensity=OFF, 
        integrationRule=SIMPSON, numIntPts=5)
    mdb.models['Model-1'].analyticalFields['AnalyticalField-2'].setValues(
        expression=k12 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-3'].setValues(
        expression=k13 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-4'].setValues(
        expression=k14 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-5'].setValues(
        expression=k15 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-6'].setValues(
        expression=k16 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-7'].setValues(
        expression=k17 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-8'].setValues(
        expression=k18 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-9'].setValues(
        expression=k19 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-10'].setValues(
        expression=k20 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-11'].setValues(
        expression=k21 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-12'].setValues(
        expression=k22 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-13'].setValues(
        expression=k23 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-14'].setValues(
        expression=k24 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-15'].setValues(
        expression=k25 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-16'].setValues(
        expression=k26 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-17'].setValues(
        expression=k27 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-18'].setValues(
        expression=k28 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-19'].setValues(
        expression=k29 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-20'].setValues(
        expression=k30 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-21'].setValues(
        expression=k31 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-22'].setValues(
        expression=k32 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-23'].setValues(
        expression=k33 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-24'].setValues(
        expression=k34 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-25'].setValues(
        expression=k35 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-26'].setValues(
        expression=k36 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-27'].setValues(
        expression=k37 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-28'].setValues(
        expression=k38 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-29'].setValues(
        expression=k39 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-30'].setValues(
        expression=k40 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-31'].setValues(
        expression=k41 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    mdb.models['Model-1'].analyticalFields['AnalyticalField-32'].setValues(
        expression=k42 + ' * ( 1.6- sqrt ( pow (  X ,2 ) + pow (  Z ,2 ) ) / 225 + 0 *  Y ) ')
    session.viewports['Viewport: 1'].view.setValues(session.views['Iso'])
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#5 ]', ), )
    region = p.Set(faces=faces, name='Set-1')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-1', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#0 #9000 ]', ), )
    region = p.Set(faces=faces, name='Set-2')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-2', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#28000000 ]', ), )
    region = p.Set(faces=faces, name='Set-3')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-3', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#0 #280 ]', ), )
    region = p.Set(faces=faces, name='Set-4')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-4', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#a ]', ), )
    region = p.Set(faces=faces, name='Set-5')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-5', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#c00000 ]', ), )
    region = p.Set(faces=faces, name='Set-6')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-6', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#3000000 ]', ), )
    region = p.Set(faces=faces, name='Set-7')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-7', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#3000 ]', ), )
    region = p.Set(faces=faces, name='Set-8')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-8', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    session.viewports['Viewport: 1'].view.setValues(session.views['Iso'])
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#50 ]', ), )
    region = p.Set(faces=faces, name='Set-9')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-9', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#0 #140 ]', ), )
    region = p.Set(faces=faces, name='Set-10')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-10', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#80000000 #2 ]', ), )
    region = p.Set(faces=faces, name='Set-11')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-11', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#0 #4800 ]', ), )
    region = p.Set(faces=faces, name='Set-12')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-12', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#a0 ]', ), )
    region = p.Set(faces=faces, name='Set-13')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-13', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#300000 ]', ), )
    region = p.Set(faces=faces, name='Set-14')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-14', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#14000000 ]', ), )
    region = p.Set(faces=faces, name='Set-15')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-15', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#c000 ]', ), )
    region = p.Set(faces=faces, name='Set-16')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-16', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    session.viewports['Viewport: 1'].view.setValues(session.views['Iso'])
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#100 ]', ), )
    region = p.Set(faces=faces, name='Set-17')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-17', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#0 #20 ]', ), )
    region = p.Set(faces=faces, name='Set-18')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-18', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#0 #8 ]', ), )
    region = p.Set(faces=faces, name='Set-19')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-19', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#0 #400 ]', ), )
    region = p.Set(faces=faces, name='Set-20')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-20', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#200 ]', ), )
    region = p.Set(faces=faces, name='Set-21')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-21', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#80000 ]', ), )
    region = p.Set(faces=faces, name='Set-22')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-22', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#40000000 ]', ), )
    region = p.Set(faces=faces, name='Set-23')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-23', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#10000 ]', ), )
    region = p.Set(faces=faces, name='Set-24')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-24', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    session.viewports['Viewport: 1'].view.setValues(session.views['Iso'])
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#400 ]', ), )
    region = p.Set(faces=faces, name='Set-25')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-25', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#0 #10 ]', ), )
    region = p.Set(faces=faces, name='Set-26')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-26', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#0 #2000 ]', ), )
    region = p.Set(faces=faces, name='Set-27')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-27', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#0 #4 ]', ), )
    region = p.Set(faces=faces, name='Set-28')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-28', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#800 ]', ), )
    region = p.Set(faces=faces, name='Set-29')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-29', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#40000 ]', ), )
    region = p.Set(faces=faces, name='Set-30')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-30', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#0 #1 ]', ), )
    region = p.Set(faces=faces, name='Set-31')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-31', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#20000 ]', ), )
    region = p.Set(faces=faces, name='Set-32')
    p = mdb.models['Model-1'].parts[partname]
    p.SectionAssignment(region=region, sectionName='Section-32', offset=0.0, 
        offsetType=MIDDLE_SURFACE, offsetField='', 
        thicknessAssignment=FROM_SECTION)
    a = mdb.models['Model-1'].rootAssembly
    session.viewports['Viewport: 1'].setValues(displayedObject=a)
    session.viewports['Viewport: 1'].assemblyDisplay.setValues(
        optimizationTasks=OFF, geometricRestrictions=OFF, stopConditions=OFF)
    a = mdb.models['Model-1'].rootAssembly
    a.DatumCsysByDefault(CARTESIAN)
    p = mdb.models['Model-1'].parts[partname]
    a.Instance(name=partname+'-1', part=p, dependent=ON)
    session.viewports['Viewport: 1'].assemblyDisplay.setValues(mesh=ON)
    session.viewports['Viewport: 1'].assemblyDisplay.meshOptions.setValues(
        meshTechnique=ON)
    p = mdb.models['Model-1'].parts[partname]
    session.viewports['Viewport: 1'].setValues(displayedObject=p)
    session.viewports['Viewport: 1'].partDisplay.setValues(sectionAssignments=OFF, 
        engineeringFeatures=OFF, mesh=ON)
    session.viewports['Viewport: 1'].partDisplay.meshOptions.setValues(
        meshTechnique=ON)
    elemType1 = mesh.ElemType(elemCode=S4R, elemLibrary=EXPLICIT, 
        secondOrderAccuracy=OFF, hourglassControl=DEFAULT)
    elemType2 = mesh.ElemType(elemCode=S3R, elemLibrary=EXPLICIT)
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    faces = f.getSequenceFromMask(mask=('[#ffffffff #ffff ]', ), )
    pickedRegions =(faces, )
    p.setElementType(regions=pickedRegions, elemTypes=(elemType1, elemType2))
    p = mdb.models['Model-1'].parts[partname]
    f = p.faces
    pickedRegions = f.getSequenceFromMask(mask=('[#ffffffff #ffff ]', ), )
    p.setMeshControls(regions=pickedRegions, elemShape=QUAD, technique=STRUCTURED)
    p = mdb.models['Model-1'].parts[partname]
    p.seedPart(size=k9, deviationFactor=0.1, minSizeFactor=0.1)
    p = mdb.models['Model-1'].parts[partname]
    p.generateMesh()
    a = mdb.models['Model-1'].rootAssembly
    a.regenerate()
    a = mdb.models['Model-1'].rootAssembly
    session.viewports['Viewport: 1'].setValues(displayedObject=a)
    session.viewports['Viewport: 1'].assemblyDisplay.setValues(mesh=OFF, 
        adaptiveMeshConstraints=ON)
    session.viewports['Viewport: 1'].assemblyDisplay.meshOptions.setValues(
        meshTechnique=OFF)
    mdb.models['Model-1'].ExplicitDynamicsStep(name='Step-1', previous='Initial')
    session.viewports['Viewport: 1'].assemblyDisplay.setValues(step='Step-1')
    session.viewports['Viewport: 1'].assemblyDisplay.setValues(loads=ON, bcs=ON, 
        predefinedFields=ON, connectors=ON, adaptiveMeshConstraints=OFF)
    a = mdb.models['Model-1'].rootAssembly
    e1 = a.instances[partname+'-1'].edges
    edges1 = e1.getSequenceFromMask(mask=('[#88 #24000040 #82000008 ]', ), )
    region = a.Set(edges=edges1, name='Set-1')
    mdb.models['Model-1'].EncastreBC(name='BC-1', createStepName='Step-1', 
        region=region, localCsys=None)
    mdb.models['Model-1'].TabularAmplitude(name='Amp-1', timeSpan=STEP, 
        smooth=SOLVER_DEFAULT, data=((0.0, 0.0), (1.0, 1.0)))
    a = mdb.models['Model-1'].rootAssembly
    s1 = a.instances[partname+'-1'].faces
    side1Faces1 = s1.getSequenceFromMask(mask=('[#f67feffc #7f7f ]', ), )
    s1 = a.instances[partname+'-1'].edges
    side1Edges1 = s1.getSequenceFromMask(mask=('[#0:2 #800000 ]', ), )
    region = a.Surface(side1Edges=side1Edges1, side1Faces=side1Faces1, 
        name='Surf-1')
    mdb.models['Model-1'].Pressure(name='Load-1', createStepName='Step-1', 
        region=region, distributionType=UNIFORM, field='', magnitude=0.4, 
        amplitude='Amp-1')
    del mdb.models['Model-1'].loads['Load-1']
    a = mdb.models['Model-1'].rootAssembly
    s1 = a.instances[partname+'-1'].faces
    side1Faces1 = s1.getSequenceFromMask(mask=('[#f67feffc #7f7f ]', ), )
    region = a.Surface(side1Faces=side1Faces1, name='Surf-2')
    mdb.models['Model-1'].Pressure(name='Load-1', createStepName='Step-1', 
        region=region, distributionType=UNIFORM, field='', magnitude=k10, 
        amplitude='Amp-1')


